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ABSTRACT 

We report on a comprehensive and consistent investigation into the X-ray emission from 
GX 339-4. All public observations in the 1 1 year RXTE archive were analysed. Three dif- 
ferent types of model - single powerlaw, broken powerlaw and a disc + powerlaw - were fitted 
to investigate the evolution of the disc, along with a fixed gaussian component at 6.4 keV to 
investigate any iron line in the spectrum. We show that the relative variation in flux and X-ray 
colour between the two best sampled outbursts are very similar. The decay of the disc temper- 
ature during the outburst is clearly seen in the soft state. The expected decay is Semsc oc T 4 ; 
we measure 7 l4 - 75±0 - 23 . This implies that the inner disc radius is approximately constant in the 
soft state. We also show a significant anti-correlation between the iron line significant width 
and the X-ray flux in the soft state while in the hard state the EW is independent of the flux. 
This results in hysteresis in the relation between X-ray flux and both line flux and EW. To 
compare the X-ray binary outburst to the behaviour seen in AGN, we construct a Disc Frac- 
tion Luminosity Diagram for GX 339-4, the first for an X-ray binary. The shape qualitatively 
matches that produced for AGN. Linking this with the radio emission from GX 339-4 the 
change in radio spectrum between the disc and power-law dominated states is clearly visible. 
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1 INTRODUCTION 

The galactic X-ray bina ry GX 339-4 (=V82 1 Ara) was discovered 
by the OSO-7 satellite jMarkert et"al]| 19731) . The optical counter- 
part to the companion star of the compact object has not been de- 
tected. Therefore the system has been classed as a low-mass X- 
ray binary (XR B) from uppe r limits on the l umin osity of the com- 
panion <Sha hbaz etai]l200lh . IHvnes et al.1 {2004) using the non- 
detection and orbital parameters, conclude that the companion is 
probably a sub-giant, of spectral type G or later. Studies of the X- 
ray spectral and temp oral characteristics c o nclude that it is a black 
hole X-ray binary dZdziarski et al.1 [T998I : ISunvaev & Revnivtsevl 
2000); the characteristics of the Fourier power spectra are more 
similar to other black hole X-ray binaries, rather than those 
of neutron stars. For other stud i es of th e X-ray emissi o n from 



this 



source see e. g 



Uedaetal.1 dl994l) 



Wilms et al.l J 1999b: iKong et al.1 fe oOO): Wa rdzinski et al 
" l2004b ' 



1998); 



2002); 



Zdziarski et alJfe004b : lBelloni et alj J2005h . 

Black hole X-ray binaries exhibit t wo main spectral 
states, the low-hard and the high-soft (e.g. Ivan der Klisl [l99l : 
iMcClintock & R emillard 2006). In a variety of accretion models, 
the only parameter which determines the state of an X-ray binary is 
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the mass accretion rate (e.g. lEsinet alJ 1997b . however a number of 
phenomena indicate that other parameters are also important. The 
major one is the hysteresis of the state transition - the transition 
from low-hard to high-soft occurs at highe r fluxes than the retur n 
to the low-hard from the high-soft state (e.g. lMivamoto et al ■119951) . 
In addition, the luminosity of the source wh en the first sta te transi- 
tion occurs varies from outburst to outburst ( Belloni 2006). 

As described in detail in Fend er et al ] d2004b . an average X- 
ray binary spends a large fraction of its time in the low-hard state. In 
this state radio emission indicates that the black hole is producing 
a steady radio-synchrotron emitting jet and the X-ray spectrum can 
be modelled by a hard powerlaw (photon index, V ~ 1.5). As the 
outburst progresses the luminosity of the binary increases and the 
X-ray spectrum softens. Eventually the jet disrupts and the X-ray 
spectrum is dominated by the soft disc emission. This is known 
the high-soft state. During the outburst, matter is accreted onto the 
black hole. The mass accretion rate reduces and disc is expected 
to cool, and therefore fade, as the outburst progresses. At the end 
of the outburst the spectrum hardens as the system returns to the 
low-hard state and the jet reforms. 

RXTE was launched in December 1995 and has been observ- 
ing GX 339-4 since lune 1996. The archive of data built up over 
the intervening period allows an extensive investigation into the be- 
haviour of GX 339-4 over a comparatively long baseline using the 
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same instrument. The long term variability allows comparisons to 
be made between the four clear outbursts seen in the RXTE data. 
The large number of individual observations also allow detailed 
studies of single outbursts of GX 339-4. 

In the hard state, GX 339-4 is a weak but steady radio source, 
with flux densities < 15 m Jy at centimetre w avelengths and a flat 
spectrum (for a review see lFender et al. 1999). In the low-hard state 
there is a non-linear positive correlatio n between the radio flux and 
both the hard and soft X-ray fluxes dCorbel et ai] |200Cl) . This is 
likely to result from a coupling between the corona and the jet. Ra- 
dio observations during the transition to the soft state have also seen 
ejection even ts which have been interpreted as larg e-scale relativis- 
ts radio jets dFender et al.lll997r . lGallo et alj|2004fl . 

GX 339-4 regularly varies by 4-5 orders of magnitude in X- 
ray flux. This relatively high number of outbursts, from a very low 
quiescent state to a high state, makes GX 339-4 one of the most 
interesting galactic X-ray binary sources. Combining radio with X- 
ray observations of the source during an outburst led to the devel- 
opment of a physical ex planation for the evolution of the sourc e as 
the outburst progresses dHoman et al. 2001; Fender et alj|2004f) . 

Black holes themselves are described by only two parameters, 
their mass and spiifl The observational properties of the X-ray bi- 
nary arise from the interplay of the parameters of the black hole 
(mass and spin) and those of the accretion flow (including the mass 
accretion rate through the disc). The black hole parameters are un- 
likely to change much on short timescales, and so it is the mass 
accretion rate and geometry of the accretion flow which is respon- 
sible for the large changes in the emission spectrum. 

Recent developments have been made in creating unifica- 
tion schemes between X-ray binaries and Active Galactic Nuclei 
(AGN). The radio-X-ray flux relation for GX 339-4 was gener- 
alised to a large number of black hole X-ray binaries bv lGallo et al.l 
d2003l) . A scaling across a wide range of mass scales from X- 
ray bi narie s to AGN wa s subse quently presented bv lMerloni et al.l 
d2003l) and lFalckeetaD d2004h and became known as the "funda- 
mental plane". Another plane using th e timing properties of b lack 
holes has recently be en constructed by McHard v et alj d2006h and 
IKording et all d2007l) . 

The observed variability timescales of AGN are many orders 
of magnitude longer than those in X-ray binaries. To study the evo- 
lution of an AGN outburst a population study is required. In order to 
be able to compare t he outbursts of black hole binaries with AGN 
IKording et all 12006) create a more general diagram to show the 
spectral states in AGN by comparing the disc flux and the power- 
law flux. We complement this by creating an equivalent diagram 
for this X-ray binary (see Section|7] 

We present our data reduction scheme in Sections [2] and [3] 
Our investigations in to the iron line and the Disc are discussed in 
Sections [3] and [6] with the Radio correlations presented in Section 

m 



2 DATA REDUCTION 

In order to fully investigate the existence and properties of the iron 
line and disc in GX 339-4, all observations in the RXTE archive 



1 Charge is a further parameter of a black hole, but as no macroscopic 
object in the Universe has as yet found to be significantly charged, we ignore 
this parameter in this discussion. 



were downloaded and analysec^]. This gave 2.59 Ms of PCA ex- 
posure over an 1 1 year baseline for this study, containing the two 
well-sampled outbursts, as well as two other outbursts with fewer 
public RXTE observations. 

All the data were reprocessed so all observations had the same 
version of the data reduction scripts applied. The HEXTE data were 
analysed so that the high energy spectral slope could be well con- 
strained. This helps in fitting the lower energy spectrum. We use the 
data reduction tools from HEASOFT0 version 6.3.2. The data re- 
duction and model fitting were automated so that each observation 
was treated in exactly the same way. 



2.1 PCA Reduction 

The PCA data were reduced according to the RXTE Cookboolfl 
and only a quick summary is given here. 

We used only data from PCU-2 as it is always switched on, 
and so can be used over the entire archive of data. It is also the 
best calibrated of the PCUs on RXTE. Background spectra were ob- 
tained using PCABACKEST from new filter files created using XTE- 
FILT, from which updated GTI files were also created. The spectra 
were extracted and the customary systematic error of 1 per cent was 
added to all spectra using GRPPHA. 

In order that the model fitting in XSPEC was reliable and rela- 
tively quick, we only fitted spectra which had more that 1000 back- 
ground subtracted PCA counts. In total there were 913 spectra ex- 
tracted from PCA data. Cutting those with fewer than 1000 counts 
from further analysis, removed 191 observations (the bottom 21 per 
cent), leaving 722. The excluded observations occur throughout the 
light curves of the outburst, with a concentration in the low flux pe- 
riods. The low-luminosity observations with low counts fall into 
the "stalks" of the Hardness Intensity Diagrams. As our analysis 
concentrates on the disc parameters, excluding these observations 
is unlikely to bias our conclusions. 



2.2 HEXTE Reduction 

The HEXTE data were also reduced according to the RXTE Cook- 
book. Where possible, we used both Cluster A and Cluster B data. 
Background spectra were obtained using HXTBACK. Spectra were 
extracted using the routine appropriate to the data-type (Event 
or Archive). Dead-time was then calculated using HXTDEAD. All 
spectra for a given ObsID were then summed using SUMPHA, 
and the appropriate responses and ancillary files were added in as 
header key words using GRPPHA. In total there were 701 spectra 
extracted from Cluster A data and 905 for Cluster B. In order to 
accurately determine the slope of the high energy power-law we re- 
quire HEXTE data to be present when fitting a model. There have to 
be at least 2000 background subtracted counts in one of the HEXTE 
clusters, with the other having a positive number of count^f]. As our 
analysis concentrates on the properties of the disc and any emission 
line which might be present, we also require PCA data to be present 
in the spectra which are fitted. 



2 The cut-off date for inclusion in this analysis was the 14 th February 
2008. 

3 |http : / /heasarc . gsf c . nasa . gov/lheasof t/| 

4 It^^^^^^^^^sf^^iasa^goWdoc^^^'^^ecipes/cookJjook . html| 

5 The background subtraction on some observations resulted in a negative 
number of HEXTE counts in one of the clusters. 
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3 MODEL FITTING 

The spectra were fitted in XSPEC (vl2.3.1ao). As our spectral anal- 
ysis concentrates on the disc properties of GX 339-4, we wish to 
analyse the spectra to the lowest possible energies. The energy 
boundaries corresponding to channel numbers have changed during 
the lifetime of the RXTE mission. The calibration of channel num- 
bers ^ 6 is uncertain, and we therefore ignore all PCA channels 
^ 6 (~ 3 keV), which allows a consistent lower bound to the spec- 
tra, extending them to the lowest possible energies and maintaining 
calibration. All PCA data greater than 25 keV were also ignored. 
The HEXTE data were fitted between 25 and 250 keV. To investi- 
gate the presence or absence of an iron line and also to see whether 
any disc emission was present, three types of models were fitted: 
Powerlaw (POWER), Broken Powerlaw (BKNPOWER) and Power- 
law + Disc (POWER+DISCBB). For each of these, a version includ- 
ing a Gaussian line fixed at 6.4 keV was also fitted, giving a total of 
six models which were fitted to each spectrum. Galactic absorption 
was modelled using the WABS photoel ectric absorption code, with 
a fixed A H = 0.4 x 10 22 cm -2 (e.g. lMiller et alj2004bh . Cabanac 
et al. (in prep.) discuss possible variations of the Ah during an out- 
burst of GX 339-4. There is insufficient spectral coverage at lower 
energies to allow the Ah to remain free during the spectral fitting. 
To obtain fluxes outside of the RXTE observing band, for the disc 
for example, dummy responses were created within XSPEC. 

To account for the galactic ridge emission, we add in the 
equivalent of a correction file to each spectrum. This has the great- 
est effect on the low flux observations. To model the galactic ridge 
emission, we used observations of GX 339-4 in Obs ID P91105 
between MJD 53550 and 53650 (25 observations, 19.1 ks expo- 
sure). We assume that the emission in these observations is dom- 
inated by the galactic ridge, with almost no contribution from 
GX 339-4. We fit this data using an absorbed powerlaw with a 
gaussian emission line. The best-fit values (F = 2.127, A/*r = 
1.153 x 10" 4 keV _1 cm 2 s _1 , line energy= 6.491 keV, a = 
1.729 x 10" 3 keV and Aline = 1-07 x 10" 6 keV _1 cm 2 s _1 ), 
are added into each model fit as a fixed set of components to take 
account of the galactic ridge emission. 

The high energy channels from HEXTE suffer from a high 
background leading to a low number of net counts in the low ex- 
posure time observations. This results in the high energy flux not 
being well determined. Using GRPPHA, we investigated the effect 
of binning. The channels were binned so that there was a minimum 
number of counts within each bin. To have clearly defined fluxes 
at all HEXTE energies we found that the level of binning had to be 
very high, and the lower energies were then severely over-binned. 
However, the change in the value of the parameters, their uncertain- 
ties and the reduced \ 2 changed very little, on the order of 5 per 
cent or less. We therefore do not perform any binning during the 
analysis of the datfl 

We encountered some difficulty in fitting discs reliably. The 
response of RXTE is only reliable down to ~ 3 keV, whereas the 
disc we are trying to measure has a temperature of around 1 keV. 
We set the minimum disc temperature to fceT = 0.1 keV in XSPEC 
to prevent discs from being fit at very low temperatures. It is likely 
that in cases where this occurred, the disc was being fit to take 
into account of any curvature in the powerlaw slope rather than to 
a true disc component. Before selecting the best fitting model we 
penalised disc models which had ksT < 0.4 keV. In doing this we 

6 In the combined spectra in Fig.|5]we bin the HEXTE points to improve 
the clarity of the spectra only. 




MJD.50000 



Figure 1. The RXTE-PCA lightcurve of GX 339-4 over 1 1 years. The flux 
is the absorbed flux between 3 — 10 keV. 

note that we are likely to have excluded a few disc fits which are 
going to be reliable. 

We initially select the model with the lowest reduced x 2 . 
However, if this best fitting model is a broken powerlaw or a 
disc+powerlaw model, we test whether using this more com- 
plex continuum model is a significant enough an improvement 
over a simple unbroken powerlaw, by performing an F-test. If 
the F-statistic probability V < 0.001 then we select the more 
complex continuum model. We investigated the number of discs 
deemed "significant" using different probability cut-offs. The num- 
ber of discs was for all intents and purposes constant regardless of 
whether the probability cut-off was 0.05 (234 disc detections) or 
(214 disc detections). We can therefore be certain that 
almost all disc values used in the rest of the analysis are from sig- 
nificant discs, what ever cut-off we use 

We note that we may be missing non-dominant disc compo- 
nents in, for example, the hard state. Even though we fit all ob- 
servations with a powerlaw + disc model, the lack of sensitivity 
of the RXTE PCA below ~ 3 keV makes detecting non-dominant 
discs difficult. Therefore even if no significant disc is detected in 
our analysis, there may be discs present at a very low level in the 
hard and intermediate states - absence of evidence does not mean 
evidence of absence. See Section[4]for some higher signal-to-noise 
spectra and further discussion of disc components when GX 339-4 
is not in the soft state. 

To determine whether the gaussian component at 6.4 keV is 
likely to be an accurate representation of a true iron line we fol- 
lowed the scheme outlined below. The line was deemed to be not 
well constrained when the width, a, was larger than 1 keV. In these 
cases, although the best fitting model includes a line component, 
only the base continuum model was used for further analysis. Sub- 
sequently, an F-test was used to determine if the addition of the 
line was statistically significant (V < 0.001), and the appropriate 
model then chosen for further analysis. 

We note that there have been recent doubts about the applica- 
bility and accuracy of using an F-test to det ermine the presence of a 
line and more complex continuum models. Prot assov et al.l (2002) 
state that the F-test does not adhere to the F distribution, even 
asymptotically. As we are attempting to determine the presence of 

7 The latter probability cut-off was chosen to be an extreme value. 
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disc and line components, this is of importance to the work pre- 
sented here. However, we note that the alternatives to an F-test are 
not easily implemented within the current fitting proceedure. Our 
solution is therefore as follows. A usual cut-off for a spectral fea- 
ture to be significant is V < 0.001 (99.9 per cent confidence). We 
investigated the number of lines deemed "significant" using differ- 
ent F-probability cut-offs. There is a larger variation in the num- 
bers of significant lines detected than in the equivalent investigation 
for the discs. Using a cut-off of 0.05 there are 533 observations with 
lines deemed significant, but using 0.00001 there are only 322. It 
is therefore not clear which F probability to use as the cut-off for 
a significant line detection. Therefore we include the estimation of 
the line significance using the normalisation of the gaussian com- 
ponent. 

This method uses the normalisation of the line component 
from XSPEC and its uncertainty. We use the uncertainty as an es- 
timate of the normalisation sigma, and hence treat line components 
as significant if their normalisation differs from zero by at least 
three sigma. Both of these two statistical tests have to be satisfied 
for the line to be used in any further analysis. Using an F-test prob- 
ability of V < 0.001 and significance of the line normalisation, re- 
sults in 400 line detections. We do note that some may be erroneous 
"detections" but the majority will be true features in the spectrum, 
see also Section|4]and Fig. [5]for higher signal-to-noise spectra. 

Any observation with a 3 — 10 keV flux from the best fitting 
model of less than 1 x 10 -11 erg s _1 was also discarded from 
further analysis, as were ones where the flux was not well deter- 
mined (the error on the 3 — 10 keV flux was larger than the flux 
itself). This resulted in a final list of 628 observations, correspond- 
ing to 2.261 Ms, with well fitted spectra and high enough fluxes 
and counts. We extract a range of parameters and fluxes for differ- 
ent bands and model components which are presented and analysed 
further the following sections. 

The full PCA lightcurve of GX 339-4 is shown in Fig. □ 
There are four outbursts covered in this period: Outburst 1 from 
MJD-50000 = 0-1500, Outburst 2 from 2340-2800, Outburst 3 from 
3040-3510 and Outburst 4 from 4000-. The MJD shown is for the 
midpoint of each individual observation. The best sampled out- 
bursts are the middle two, and these are the ones on which some 
of the remaining analysis concentrates on as these have the most 
observations. 

F o llowing the analysi s in e.g. lHoman et al. | l l200ll) : iBellonil 
( 2004): iFender et all |2004), we plot the X-ray observations in a 
Hardness-Intensity Diagram (HID) to show the state changes of 
GX 339-4 during its outbursts. We extract the 3 — 6, 6 — 10 and 
3 — 10 keV fluxes from the spectra. The X-ray colour was calcu- 
lated from S^iokcv / '^s-skeV and plotted against S3-iokeV in 

Fig.m 

In Fig. [2] we show the X-ray colours of GX 339-4 where the 
state transitions between the high-soft and soft-intermediate states, 
as well as the hard-intermediate and low-hard states occur. The 
transitions have been determined by changes in the ti ming prop- 
erties of the source as calculated in lBelloni et alj d2006t) using data 
corresponding to our Outburst 3 (see Section [4}. Both transitions 
towards and away from the soft state occurred between observa- 
tions such that a single X-ray colour could be used to delineate the 
state changes. 

iBelloni et alj J2005h study our Outburst 2 and the transitions 
from the hard state and into the soft state are very similar to those 
in Outburst 3 and so match those we use in this work. However, 
the transition between the hard- and soft-intermediate states occurs 
at a different X-ray colour. Also, there is no observation in soft- 
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Figure 2. TOP: The hardness-intensity diagram from all the observations. 
Powerlaw models in squares, Broken Powerlaw models in circles, Pow- 
erlaw + Disc models in triangles - i.e. all triangular points require discs. 
X-ray colours of 0.225, 0.41 and 0.85 define the soft-soft-intermediate, 
soft-intermediate-hard intermediate and hard-intermediate-hard state tran- 
sitions. BOTTOM: The HIDs using a single continuum model only powerlaw 
(PL), broken powerlaw (BPL) and disc + powerlaw (DPL), left to right. The 
axes are the same scale as the main figure. 



intermediate state on the return to the low-hard state in Outburst 2. 
The X-ray colour for the transition from the hard-intermediate to 
the soft-intermediate state in Outburst 3 is S , 6-iokev/<S , 3_6kcV < 
0.37. 

There is no clear reason why the transition should occur at 
the same X-ray colour for any two outbursts. The X-ray colour is 
determined by the temperature and photon index, and also by the 
relative strength of the two components. The transitions towards the 
soft state are very different between the two outbursts, and so we 
should not expect the same X-ray colour to work for all three tran- 
sitions. As a single X-ray colour could be determined for the tran- 
sition between the two intermediate states in Outburst 3 we show 
this on the appropriate figures, but do not show an equivalent no 
figures for Outburst 2. 

Initial investigations into the HID for the most recent outburst 
(Outburst 4) whose transition into the soft state occurred at a similar 
flux to Outburst 2, shows transitions between states calculated from 
X-ray timing, at similar X-ray colours to Outburst 2 (Del Santo et 
al., in prep.). 

From now on we concentrate main ly on the soft and the hard 
states as these are the best sampled (e.g. lBellomll2006l) . For the re- 
mainder of this work we define the soft state when the X-ray colour, 
S6-10/S3-6 < 0.22, and the hard state when Se- 10 /S^-e > 0.87 
unless stated otherwise (see Fig.[2l, and the intermediate state is not 
split in any discussions. 

We also show the HIDs obtained for GX 339-4 using each of 
the thee continuum models individually in Fig. [2] The basic shape 
of the HID is recovered in all three cases, with the most similar be- 
ing the one for the disc + powerlaw. There is a clear "pile-up" of 
observations with an X-ray colour around 0.2, which is close to our 
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adopted transition into the soft state. This is mainly the result of 
our restriction on the slope of the powerlaw to have V < 4.0. This 
limit was set to aid the fitting of a disc. As the PCA is only sensi- 
tive down to ~ 3 keV the curvature of the disc component can be 
difficult to fit when the disc is not dominating the spectrum. There- 
fore a fit might result in an extremely steep powerlaw rather than a 
disc. Also, if there are insufficient HEXTE counts then the power- 
law slope cannot be well determined from the hard X-rays, which 
would allow an nonphysical steep powerlaw fit with the disc in the 
soft state. Powerlaws of Y > 4.0 are thought to be unphysical. The 
best fitting models just harder of this X-ray colour are ones with a 
disc component, and observations in the intermediate states which 
have broken powerlaw models as their best fit have well determined 
slopes not pegged at 4.0. Therefore we do not believe that this re- 
striction is influencing our results (see Section[4]for disc models in 
the intermediate state). There is little difference between the pow- 
erlaw and broken powerlaw models as the break usually occurs at 
around 10 keV and so has little effect on the X-ray colour adopted 
here. 



3.1 Comptonisation Models 

We note that the powerlaw models we use here have limited physi- 
cal meaning. We therefore performed a quick investigation into fit- 
ting the COMPTT model combined with a DISKBB and gaussian line 
component to the data. Using this combination of models would re- 
sult in a disc component being fitted whatever the state of the source 
and so would be a method for investigating the disc parameters out- 
side of the soft state. Selecting the four hard state spectra with the 
greatest number of PCA counts resulted in a good fit with well de- 
fined values for the r and kT parameters, which were reasonably 
similar between different observations (~ 2.5 and ~ 25 keV re- 
spectively). 

However, turning to the soft state, again using four spectra 
with the largest number of PCA counts, we found that the lack of 
high energy counts from HEXTE meant that we could not determine 
both r and kT for every observation. Fixing kT — 50 keV, the 
average value from those observations where both parameters could 
be estimated, allowed a value of r ~ 0.5 to be estimated. 

Using a single analysis scheme and script for all the observa- 
tions is vital to remove any systematic differences between the anal- 
ysis of different observations. As the values of both r and kT varied 
between the soft and hard states, and the fact that we were not able 
to determine independent values for both parameters in the observa- 
tions of the soft state with the highest PCA counts, we did not inves- 
tigate any further. The observations investigated with COMPTT had 
between 0.7 and 1 1 million PCA and 60 and 350 thousand HEXTE 
background subtracted counts, whereas our cut-off for the number 
of counts in an observation is 1000 PCA and 2000 HEXTE back- 
ground subtracted counts. We do not believe that many well defined 
parameters will be extracted from fitting COMPTT models to these 
shorter, low counts observations. 



3.2 Best Fitting Models 

The description of the method used to select the best fitting model 
is outlined in Section [3] Fig. [2] shows which was the best fitting 
model through the HID. Models without the disc component are 
the best fitting for the majority of the hard and intermediate state 
observations. The soft state observations are exclusively best fit by 
models which include a disc component. 



However, some of the observations in the hard state have best 
fitting models which include a disc component. It is expected that 
the disc emission should be very low compared to the power-law 
emission when GX 339-4 is in the hard state. As a result it is sur- 
prising that for some hard state observations the best fit model in- 
cludes a disc component. It is likely that, although these models 
have the lowest x 2 from the ones fitted, their disc parameters are 
not accurate. These disc fluxes are one to two orders of magnitude 
below those from observations in the soft state. The errors on the 
disc temperature from these fits are larger than those obtained from 
observations in the soft state. A combination of a broken power- 
law and the galactic absorption may have combined to emulate the 
emission from a disc sufficiently that a disc model is the best fitting. 

Other X-ray missions (e.g. Swift, XMM-Newton) have a re- 
sponse at low energies such that, if a disc were present in the hard 
state, then they would be able to detect it in a suitable length ob- 
servation. In GX 339-4 Cab anac et al. (in prep .) detect discs out- 
side of the soft state, as do iMiller et al However, in the 
clos er source XTE Jl 1 1 8+480 discs have been seen in the hard 
state lChatv eTal] d2003l) : lMcClintock et all feOOll) . 

The distributions of the reduced for the three best fitting 
continuum models are shown in Fig. [5] The range in \ 2 f° r me 
three continuum models is on the whole between 0.8 — » 1.5 with 
some outliers. Most outliers are in the simple power-law model fits 
as this is the default model, but they are few in number. 



4 THE TWO MAJOR OUTBURSTS 

We give a short summary of the major features of the two main 
outbursts of GX 339-4. Most of the observations in the rest of 
this study (424/634) come from these two outbursts. The detailed 
lightcurves, along with the X-ray colours and HIDs are shown in 
Fig. E] Outburst 2 reaches its peak flux on the initial rise. There 
is a small decay as GX 339-4 evolves into the soft-state, and then 
the flux rises again to the same maximum flux. After this there is a 
gr adual decay in flux. A detailed analysis of this outburst is given 
in lBelloni etail fc005h . 

In Outburst 3 the initial peak in the flux is almost an order of 
magnitude below the eventual maximum of the outburst. GX 339-4 
remains in the hard state for over 100 days after this first peak, with 
the flux falling by a factor of two, then GX 339-4 rapidly moves 
into the soft-state. The flux continues to increase once in the soft 
state, and spectrum hardens fractionally. After the break in the light 
curve, by which time GX 339-4 is back in the soft state, the spec- 
trum hardens before the flux falls and GX 339-4 rejoins the hard 
branch. This motion can be traced from the change in colour of the 
points in the HIDs shown in Fig. [4] A detailed ana lysis of the early 
part of this outburst is given in Bello niet al]j2006l) . 

It is clear that Outbursts 2 and 3 reach different peak fluxes 
on the hard branch (see Fig. [4}. Even though there are only a few 
observations, Outburst 4 also appears to peak at similar fluxes to 
that of Outburst 2 (see Fig.[T] Outburst 1 is insufficiently well sam- 
pled to be able to comment on its peak flux. On the return to the 
ha rd branch, a l l of th e outbursts are at very similar fluxes, as noted 
bv lMaccarond d2003l) for black hole X-ray binaries in general. It is 
possible that Outburst 3 is a factor 2 higher, but the offset is small 
compared to offsets in the flux of the intermediate states for each 
outburst. 

The variation of the X-ray colour during Outbursts 2 and 3 
is similar even though they occur two years apart. The change in 
X-ray colour with the outburst is shown again in Fig. [6] rescaled 
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Figure 3. The distributions of the reduced \ 2 f° r the three different continuum models selected as best fitting using the method presented in the text. The inset 
in the powerlaw model is an enlargement of the peak of the distribution. 
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Figure 4. The two best covered outbursts of GX339-4. TOP: Detailed light curves and X-ray colour plots; BOTTOM: HIDs of the two main outbursts. For 
clarity the errorbars are only shown in the flux and X-ray colour light curves rather than on the HID. On the HIDs we show the line tracking GX 339-4 from 
observation to observation. The colour scale is the date since the beginning of the outburst, blue is early and red is late. The boxes in the bottom diagrams show 
the selected outbursts for the spectral stacking described in Section^and shown in Fig|3] The number of the box corresponds to the number of the spectrum 
shown in Fig[5] The dotted vertical lines show the X-ray colours adopted for the state transitions outlined in Section[3] 



so that the beginning and end of the transition to and from the soft 
state are aligned. The "stretching factor", the relative lengths of the 
two outbursts, is around 4 : 3 for Outburst 2 to Outburst 3. The 
two curves are remarkably similar - the rise in flux and associated 
hardening just before halfway through an outburst occurs in both. 



There is also a suggestion that the temporary softening of the spec- 
trum as GX 339-4 hardens again, which is clearly seen in Outburst 
2, is also present in Outburst 3. This dip is also seen in the tail end 
of Outburst 4 when matched to the previous outbursts. 

We investigated whether disc components could be present in 
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Figure 5. The spectra and model fits to the summed spectra. The fits are for those models which have been used in the automated fitting of the individual 
spectra, except for the Low Intermediate state, where the DISCBB+POWERLAW model is shown. The more complicated models detailed in Tableware not 
shown here. The locations in the HID from where the spectra were taken are shown in Fig. [4] and are also indicated by the numbers in the top right corner of 
each spectrum. The data have been binned for the clarity of the plots alone (no binning was done prior to fitting) so that the significance of each data point is 
at least 3<r. The effect of the binning is most obvious at the high energy end of the HEXTE spectra as well as the shorter exposure time data. 
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observations where individual spectra were unable to indicate any 
significant disc component in the intermediate and hard states. We 
selected six areas of similar X-ray colour and flux from the HIDs 
of these two well observed outbursts corresponding to low hard, 
high hard, very high soft, high soft, low soft and low intermediate 
states (see boxes in Fig.©. These areas were selected to cover a 
range in states, but also where there were a concentration of obser- 
vations at a similar epoch. We sum the PCA and HEXTE spectra 
and associated backgrounds and refit the same sets of models as 
before in XSPEC (see Fig[5]for the spectra and the Appendix for the 
numerical results of the fitting). 

The soft states are, unsurprisingly, well fit by a DISCBB + 
POWERLAW model. The best fit model to the intermediate state 
spectrum is also a DISCBB + POWERLAW model. We note that this 
is different to the best fit model of some of the individual spec- 
tra which make up this summed spectrum (compare Figs. [2] and 
|4j. The improvement in \ 2 is small, but significant as the degrees 
of freedom remain the same. We therefore conclude that the im- 
proved signal-to-noise of the summed spectrum allows an accurate 
fitting of the disc and so this model wins out over a purely broken 
powerlaw model. The low intermediate summed spectrum was con- 
structed from observations of GX 339-4 as the outburst decayed, 
and as such a faint disc would be expected in this state. The disc 
temperature does, however, appear to rise relative to the soft state 
rather than fall as would be expected (see Section [6] for more de- 
tails). 

Both summed spectra from the hard state are poorly fit by the 
simple absorbed broken powerlaw model. The high flux of GX 339- 
4 during the observations or the long equivalent exposure of the 
summed spectra has resulted in very high signal-to-noise. Conse- 
quently the simple model is no longer a good fit to the data. An 
improved, but by no means good, fit to the data is achieved by in- 
troducing another break in the powerlaw slope as there is indica- 
tion for this in the HEXTE data. We also added in a DISCBB model 
just in case a disc were present and detectable in these high quality 
spectra. There was no improvement to the fit in the low hard state 
(which is taken from observations on the rise of the outburst). In the 
high hard state, however, adding a disc improved the fit (see notes 
for Table H. 

The summed spectra do show disc components when they are 
not significantly detected by all the individual observations. This 
study can therefore not shed much light on the rise or fall of the 
disc in the hard state of GX 339-4. However, in long exposure ob- 
servations in the high hard state, the fitted disc is at higher disc 
tem peratures t han w hat is observed in the soft state. 

lYu et al.l d2007h show that the peak hard state flux of an out- 
burst depends on the length of time since the previous outburst of 
GX 339-4. Using the detailed light- and X-ray colour-curves from 
the four outbursts observed by RXTE. The duration of Outbursts 2 
and 3 (their Outbursts 6 and 7) appear to correlate with the length of 
time since the previous outburst (see Table[]}. However Outburst 4 
lasts longer than would be expected from the previous two and also 
reaches the same peak flux as Outburst 2 without such a long period 
of quiescence preceding it. Under the assumption that the duration 
of the outburst is linked to the "waiting time" by a power-law re- 
lationship, then Outburst 4 would be ex pected to last around 280 
days (Outburst 4 was not in the study of lYu et aL 2007). However, 
the end of the outburst is well monitored by RXTE and there are 
sufficient observations at the beginning to rule out a much earlier 
start (a mis-measurement of 100 days would be required to match 
the relation obtained for the other outbursts) or a lower peak flux. 

The length of the outburst, however, would depend on how 
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Figure 6. The variation in X-ray colour as the outbursts progress for the 
two major outbursts observed by RXTE. The date scale as been adjusted 
to align the beginning and end of the outbursts. The colour scale is that of 
the 3 - 10 kcV flux of GX 339-4, white for high flux and black for low 
flux. The colour scales have not been adjusted to match between the two 
plots. The dotted lines show the X-ray colours adopted for the different 
state transitions outlined in Section[3] 



Table 1. Outburst Durations 



OB 


Start 


Stop 


^outburst 


^prior 


Peak Flux 


C 




MJD-50000 


days 


days 


( erg s _1 cm -2 ) 




1 


800 


1240 


440 




2.7timesl0- 9 




2 


2400 


2730 


330 


1160 


9.4 x 10" 9 


123 


3 


3225 


3475 


250 


495 


5.0 x 10" 9 


99 


4 


4130 


4240 


110 


655 


9.7 x 10~ 9 


68 



The start and stop times for the outbursts are taken as the departure and 
return to the hard state (X-ray colour =1). The separation is then the time 
between the return to the hard state for one outburst on its decline and the 
departure from the hard state on the subsequent outburst on its rise. 

quickly the material in the disc was drained during the outburst as 
well as the amount of material which has built up in the disc since 
the last outburst . Under the assumption of a constant deposition 
rate of material into the disc and that the disc is completely emptied 
during each outburst, then the interval is a measure of the amount of 
material which could contribute to the outburst. In reality the disc 
is unlikely to be totally emptied during each outburst. 

As is clear from Table [T] the length of the outburst, t ou tburst, 
does not correlate well with the interval between outbursts, t pl i OI . 
However, the duration of an outburst depends not only on the 
amount of material in the disc, but also on the rate at which it is 
depleted. The peak X-ray flux is likely to be a good measure of 
the depletion rate within the outburst. If we correct the interval be- 
tween outbursts by dividing by the X-ray flux, we obtain values 
which do correlate with the length of the outburst, C. We note that 
we only have four outbursts and hence three estimate on the relative 
magnitudes of the durations of the outbursts to use in this study. 



5 IRON LINE 

An F~test in combination with the line normalisation (see Section 
[3] for more details) were used to test for the presence of an emis- 
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sion line (assumed to be the iron K a fluorescence line) at 6.4 keV. 
However out of 628 spectral fits, 400 require a and these oc- 
cur all over the HID, in all states. 

We did investigate the effect of allowing the line energy to be 
free, between 6 — 7 keV. However this caused some erroneous fit- 
ting, mainly in cases where there is some slight curvature to the 
low energy powerlaw. The curvature is a broad feature usually at 
lower energies than 6 keV. In trying to fit this feature, the gaussian 
component became very broad and "pegged" at the lower energy 
limit with a large uncertainty in the line energy, even though the 
residuals showed clear excesses in the 6 — 7 keV range. Selecting 
these models with a significant "line" detection as the best fitting 
would not give information about any line present as this has not 
been fitted by the gaussian. Therefore, any line component, even if 
deemed significant by the F-test which had an energy uncertainty 
of > 1 keV was taken not to be representing a line, but another con- 
tinuum component, and so the line "detection" was discarded from 
further analysis. In some cases a line was well and sensibly fit at en- 
ergies different from 6.4 keV, however these were few. However, 
this method discards observations where a line may be present but 
XSPEC has preferentially fitted the curvature in the powerlaw in- 
stead. We therefore present the results on the analysis with a gaus- 
sian component fixed at 6.4 keV in this section. We note that using 
the most stringent cut-off in the F-test when selecting lines when 
their energy was free did show a variation in the energy across the 
HID. The lowest energy was at the top of the low-hard state and 
the highest energy was in the intermediate stat^f] at the end of an 
outburst, during the return to quiescence. 

We extract the line flux over a 0.1 — 10 keV range. The line is 
strongest at the highest fluxes in the hard state, see Fig. [7] whereas 
the EW is largest at the lowest flux soft states. The line width is 
also largest in the soft state (Fig.[7). Fitting the line flux against the 
3 — 10 keV total flux shows a relationship for all of the data of 



So. 



oc S. 



3-10 keV i 



however there is a large scatter of around the best fitting line (see 
Fig H}. The Spearman Rank correlation coefficient is 0.860. We 
also calculate the Kendall's Tau correlation coefficient, 0.678. The 
variance of the lattei0 is 0.0011 and so the significance of the 
correlation is 20.2<j. As the continuum increases, then the line 
does in an almost linear relationship. However, the X-ray colour of 
GX 339-4 for each observation shows that the hard and soft states 
have different behaviours. 

We split the data into extremely soft (X-ray colour < 0.22) 
and hard (X-ray colour > 0.87) states and refit the correlation. The 
hard and soft state correlations for both the 7 — 20 keV and 3 — 
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8 This is using both the .F-test and the line normalisation test to check that 
the gaussian component is significant. 

9 Given the problems with fitting a line whose energy was free, we do not 
split the intermediate state for this comment. 

10 Kendall's Tau is more non-parametric than the Spearman Rank corre- 
lation as it uses only the relative ordering of the ranks. The variance is 
calculated as Var(r) = g^^— i) 1 w ' m tne significance, m °". as then 
r/Var(r). 



Figure 7. The hardness-intensity diagram from all the observations where 
a significant line is detected. TOP The colour scale is from the line flux. The 
decrease in line flux during the course of the outburst is clearest at the low 
soft states. MIDDLE The colour scale is for the equivalent width of the iron 
line (logarithmic, eV). In some cases the calculated equivalent width was 
greater than 1 keV and these have been truncated in the colour scale to aid 
clarity. BOTTOM The colour scale is for the line width in keV. 
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Figure 8. The flux of the 6.4 keV iron line as a function of TOP the 3 — 
10 keV total flux and BOTTOM the 7-20 keV total flux. The colour scale 
on the plot is the X-ray colour of the observation. The best fit relation is 
to all the data is shown by the dotted line. The best fit to the soft and hard 
states are shown by the black solid and dashed lines respectively. The insets 
show the approximate motion of GX 339-4 through the diagrams during an 
outburst. 



10 keV fluxes are then 
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respectively (see Fig[5}. We also show the cor relations of t he 7 — 
20keV total flux with the line flux, following Ikossi et al.l J2005t) . 
as only photons with energies > 7keV are be able to ionise iron. 
We do not, however, fit all observations with a single line as there 
are clearly two separate relations in the 7 — 20 keV flux in Fig. [8] 
In both flux bands, the slope of the relation of the line flux 
is different between the soft and hard states. In the 7 — 20keV 
flux band, the powerlaw is the dominant spectral component being 
probed. The powerlaw flux falls away during the soft state and, 
noting that the lineflux remains approximately constant across all 
X-ray colours in Fig. [7j the soft state points have been displaced 
to the left. However using the lower energy 3—10 keV band, the 



disc is the dominant spectral component in the soft state, but the 
powerlaw takes it's place in the hard state. Therefore there is no 
sideways displacement of the soft state points relative to the hard 
state. 

The difference between the behaviour of the line fluxes in the 
hard and the soft state could arise from changes in the ionization 
parameter as well as from differences in the spectral slope of the in- 
cident radiation. The X-ray spectrum above 7.1 keV of a soft state 
object is steeper than that of a hard state object. Thus, there are 
more lower energy photons available to ionize the iron in the soft 
state compared to the hard state. As the cross-section of the ioniza- 
tion process decreases with increasing photon energies, a soft state 
is therefore more efficient to produce the iron line than a hard state 
for a given 7-20 keV flux. Additionally to this effect a soft state may 
have a smaller ionization parameter, which wo uld also increase the 
obser ved line flux (for more details see e.g. iRevnolds & Nowakl 
120031) . 

5.1 X-ray Baldwin Effect 

Having detected lines in most of observations at different states 
during out an outburst we i nvestigated whe ther there is an X-ray 
Baldwin effect in GX 339-4. iBaldwir] dl977h found that the equiv- 
alent width of the Civ line decreased with increasing UV lumi- 
nosity in a sample of AGN. An X-ray Ba ldwin effect has been 
repor t ed in AGN by a number of authors dlwasawa & Taniguchil 
19931: iNandra et al l 1 19971 ; IPage et alj 12004 ; iBianchi et all 120071 : 
Mattson et al.ll2007D . They find that the equivalent width depends 
on the luminosity in the following range 

EW oc I/ 3 , 

where -0.20 < < -0.14. 

The equivalent width was calculated during the course of the 
fitting, and is shown on the HID in Fig. [7] The HID clearly shows 
that, on average, the equivalent width is highest in the softest and 
lowest flux states. Fig.[9]shows the equivalent width of the iron line 
against the flux of the observation, with the X-ray colour of the 
spectrum as a colour scale. The two states (hard and soft) appear to 
have different correlations between the flux and equivalent width. 
Extracting those observations which have an X-ray colour of < 
0.22 (the extreme soft state), the best fit to the anti-correlation is, 

FW rv <?-0.61±0.09 

"3-10 keV 
BW ~ Q-O.55±0.05 

E W oc 6 7 _ 20koV , 

with Spearman Rank correlation coefficients of —0.840 and 
-0.430 respectively. The Kendall Tau are -0.644 and -0.308, 
with significances of 10.6cr and 5.07(7 respectively. We also take 
the hard state observations (those with X-ray colours > 0.87) we 
fit the observations with 

TPUf ~ q0.32±0.07 
EW OC »!53_ 10keV 

PW ~ q0.33±0.02 
E W OC ^7_ 20 keV . 

with Spearman Rank correlation coefficients of 0.262 and 0.237 
respectively. The best fitting line has a positive index. However 
this is the result of the combination of the large scatter in the 
majority of the hard state points and the small cluster at high 
fluxes, separate from the rest. If we ignore those points with S > 
2 x 10" 9 ergs _1 cm~ 2 and EW > 400 eV, and fit the bulk 



of the hard state observations, we find EW oc S„ 



and 



EW oc Sf^ofev' 04 > shown as the dotted line in Fi g- l2l These 
are consistent with no change in the EW with flux. No errors in the 
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Figure 9. The Equivalent Width of the Iron line against the flux during the 
outbursts, TOP: 7 - 20keV flux, BOTTOM: 3 - lOkeV flux. The colour 
scale shows the X-ray colour of the spectrum for each observation, black 
being soft and white being hard. The solid line shows the best fit to the soft 
state, the dashed line the best fit to the hard state, and the dotted line to the 
low-flux hard states. The arrows indicate the motion of GX 339-4 through 
the diagram during an outburst. 



7 — 20keV fluxes could be obtained from XSPEC and so the un- 
certainties in the slopes for the 7 — 20 keV fluxes only include the 
errors in the EWs. The uncertainties in the slope of the 3 — 10 keV 
fits include the errors from the EWs and the fluxes. 

If the lineflux decreases at the same rate as the continuum then 
the measurable equivalent width of the line decreases as the flux 
decreases. This effect would therefore lead to a lower cut-off for 
the minimum measurable equivalent width with flux with a negative 
slope. This would result in a small slope in the data points, whereas 
it should be flat. We therefore do not believe that there is an X-ray 
Baldwin effect in GX 339-4 in the hard state. 

The motion of GX 339-4 through the flux-equivalent width 
plane is that at the beginning of the outburst the flux increases at 
constant equivalent width. At the peak flux, when the source moves 
into the soft state, the equivalent widths start to increase as the flux 
decreases. When the source returns to the hard state, the equiva- 
lent widths reduce suddenly and then become constant as the flux 
falls (see arrows in Fig. [9}. This is another example of hysteretical 
behaviour from GX 339-4. 

There is an X-ray Baldwin effect in the 6.4 keV iron line in the 



soft state of GX 339-4. The slope of the relation, however, i s much 
steeper than that observed in AGN. Ilwasawa & Tani guchi i 19931) 
studied Seyfert Is and quasars, whereas in a sam ple of 53 AGN 
from the XMM-Newton archive . I Page et al" I J2004h used 14 radio 
loud and 39 radio quiet sources. Bia nchi et ajj J2007t) excluded ra- 
dio loud object £3 from their complete study of radio quiet Type 1 
AGN in the XMM-Newton archive. These studies all show an X-ray 
Baldwin effect between a variety of sources, rather than an intrinsic 
effect within one source. Most of the AGN in studies of the X-ray 
Baldwin effect have been radio quiet AGN, analogous to the soft 
state of black hole X-ray binaries. Scaling the effect described here 
for GX 339-4 up to AGN would result in an effect which occurs on 
very short timescales. 



5.2 Links to the Powerlaw Slope 

In a re cent study of 12 Seyfert 1 and 1.2 galaxies. iMattson et al.l 
( 2007) fitted RXTE data with a PEXRAV model, which simulates 
an exponentially cut-off power-law reflected by neutral matter, and 
a Gaussian line. They find a positive correlation between the line 
equivalent width and the underlying power law slope for T < 2. 
For r > 2 they found an anti correlation, with a peak equivalent 
width of ~ 250 eV, see their Figs, la and 3. We plot the powerlaw 
slope (below the break in the case of a broken powerlaw) against the 
equivalent width for GX 339-4 in Fig.|To] We note that superficially 
the shape of the distribution is very similar, however there is a larger 
scatter in the shape of the distribution. The colour scale corresponds 
to the X-ray colour of the observation. The correlation corresponds 
to the hard state whereas the anti-correlation corresponds to the soft 
state. 

In GX 339-4 the hard state, which usually associated with the 
pr esence of a steady, m ildly relativistic jet, is found at the low Ts. 
In IMattson et alj J2007t) the low T end of relation is dominated 
by points from 3C 273. The positive correlation is composed of 
observations of radio loud Seyfert 1 galaxies, the peak by radio 
quiet Seyfert Is, and radio quiet Seyfert 1.2s comprise the anti- 
correlation part of the relation. In GX 339-4, the soft, disc domi- 
nated state is mostly at the peak of the relation with the intermedi- 
ate state at the high F end. 

The model o f an X-ray binary outburst presented in 
iFender et ail J2004I) indicates that the jet is quenched as the disc 
becomes dominant in the soft state. As the jet is the source for the 
radio emission, when this quenches, the X-ray binary becomes ra- 
dio quiet, in the soft state, as shown in Fig. [16] We therefore would 
expect to see observations of the soft, radio quiet state of GX 339- 
4 on the negative correlation and the hard, radio loud state on the 
posi tive correlation from this model. 

iGeorge & Fabianl Jl99lh simulated the spectrum expected 
from an X-ray source which illuminates a slab and showed that 
it should include an iron line and a "Compton hump." As the spec- 
trum softens (r increases) there are fewer photons with enough en- 
ergy to photo-ionise iron, and so the iron line EW decreases. This 
is observed in the so ft state of GX 339-4 as well as in the sample in 
IMattson etail J2007l). 

IMattson et alj 120071) suggest that for F to increase as the jet 
dominance decreases the jet would have a hard X-ray component 



11 Quasars with log(i?) > 1 and Seyferts with log(ij) > 2.4 
and log(_Rx > —2.7 55) where Ft is the radio-loudness parameter 
fromlStocke et alj 1 19921) a nd Rx is the X-ray radio loudness parameter 
iTerashima & Wilsonl feOOlhjBianchi et alj2007l) 
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Figure 10. TOP: The equivalent width against the power law slope. The 
colour scale shows the X-ray colour of the spectrum for each observation, 
black being soft and white being hard. The plot has been truncated at T = 
1.4 and T = 3.5, and T outside of this range is unlikely to be an accurate 
description of the data. BOTTOM: The disc fraction of GX 339-4 against the 
equivalent width with the colour scale representing the X-ray colour. See 
Section|7]for more details about the disc fraction. 



associated with it. On the other side of the relation, observations of 
MCG-6-30-15 hav e shown that the refl ection component remains 
relatively constant jMiniutti et al.ll2007h . so when T increases the 
EW decreases. 

As the relations between the power-law slope and the equiva- 
lent width of the iron line are similar for GX 339-4 and the sample 
of AGN it is reasonable to assume that the processes at the centre 
of an AGN and an XRB are similar. As XRBs have radio (jet) dom- 
inated quiescent periods and disc dominated outbursts, then this 
similarity in the relations points towards to AGN also having out- 
burst phases. As the black hole masses in AGN are ~ 10 5 — 10 s 
times larger, their outburst phases are much longer. However, this 
result indicates that the path of an AGN through an HID would be 
similar to XRBs, though their direction is currently unknown. The 
AGN HID would have to be constructed so that ultra-violet emis- 
sion from the disc would be included. This is further discussed in 
Section|7] 

As the hard, intermediate and soft state points appear in dif- 
ferent locations in the diagram, we reform the relation using the 
disc fraction (see Section [7] and Fig, 1 10b. The X-ray colour is in- 



directly related to the disc temperature and luminosity. The disc 
fraction, however, is a more physical quantity describing the rela- 
tive luminosities of the disc and the powerlaw components of the 
spectrum. The hard and intermediate state appear with no disc frac- 
tion as the best fitting model contains no disc component. The soft 
states occur over a range of disc fractions and equivalent widths. 
We split the observations into two sets; those with a disc (disc frac- 
tion < 0.5) and those without (disc fraction ^ 0.5) and perform a 
t-test to indicate whether the means are significantly different and 
also a K-S test to see if the two distributions of equivalent width 
are drawn from the same population. The probabilities for both 
tests show that the two data-sets have the same mean/were drawn 
from the same population are very low (Pks = 1-29 x 10 -25 and 
Pt = 1.3 x 10 -46 respectively). We therefore conclude that, al- 
though the separation of the two populations is small, the behaviour 
of the equivalent width of the line is different when GX 339-4 is in 
outburst to when it is in quiescence/the hard state. However the dif- 
ference in behaviour is less clear in this presentation of the variation 
of the equivalent width with the state of GX 339-4. 



6 DISC TEMPERATURE 

Selecting those observations where a disc is present we extract the 
disc temperature. The variation in disc temperature through the 
HID is shown in Fig. QT] The highest disc temperatures are ob- 
served at the highest fluxes in the soft state, and the lowest temper- 
atures at the lowest fluxes. The variation through the soft state is 
smooth, even though we show data from all four outbursts. The 
observations with an X-ray colour > 0.4 show a disc tempera- 
ture much higher than those in the soft state (though the colour 
scale is truncated at 1 keV). In fact they appear to indicate that the 
disc temperature is rising as the source is fading down the hard 
branch. As stated in Section l3~2l it is unlikely that these disc "de- 
tections" are true. Without spectral information below 3keV we 
are limited to how accurately we can determine the presence or 
absence of a disc. If there is any curvature in the spectrum, this 
will be better fit by a higher temperature disc combined with a 
powerlaw rather than a simple powerlaw. However, the curvature 
may not be the result of a disc component, and as we have on pur- 
pose restricted the choice of models to a limited range, the disc 
turns out to be the statistical best fit. The temperature of these 
disc fits is much higher than would be expected for a decaying 
disc (> 1 keV). If these are true discs, then the temperature at the 
end of an outburst doubles over 10-20 days compared to the grad- 
ual ~250 day decay in the soft state. Also, as the fluxes of these 
disc components is 10~ n — 10~ 10 erg s _1 cm~ 2 rather than the 
10" 9 - 10~ 8 erg s cm for the soft-state disc detections, we 
conclude that the DISKBB model is not fitting the true disc compo- 
nent. We do not include any of these points in any analysis in this 
section. 

The spectrum of an optically thick, geometrically thin accre- 
tion disc is the result of a sum of black-body spectra, one for each 
radius - a so-called multi-coloured disc model. The maximum tem- 
perature for such a disc should occur close to the inner edge of the 
accretion disc, R[ n . In a similar way to the Stefan-Boltzmann law, 
the resulting disc flux should follow: 



Sbis 



oc R? n T 4 , 



mimicking the emission from a black body whose radius is that of 
the inner edge of the disc and whose temperature is that of the disc 
at the innermost radius. In the first instance we assume a constant 
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Figure 11. The hardness-intensity diagram from all the observations where 
a disc is detected. We do not consider the apparent detections of the disc in 
the hard state as true detections. The colour scale is for the disc temperature 
in keV, and has been truncated at 1 keV to emphasise the variation in disc 
temperature during the soft state. 

inner radius for the disc. We initially assumed that in the soft state, 
the flux from the disc would dominate the total flux to such an 
extent that this relation would hold for the total flux received from 
GX 339-4 and hence expected that S3-10 k c v oc T 4 . 

Fitting the correlation of the total flux from GX 339-4 with the 
disc temperature (for disc temperatures < 1 keV), we obtain 

O rri9.44i0.12 

03-10 keV <X 1 , 

which is steeper than what is expected and is shown by the dashed 
line in Fig [12] TOP. Our assumption that the flux is dominated by 
the disc component is unlikely to be true. We see extra flux above 
what would be expected from the disc alone at the lower flux obser- 
vations. As a result we need to use the unabsorbed disc flux. There 
are two methods of obtaining the disc flux from the fit; direct from 
XSPEC by setting the normalisations of all other emission compo- 
nents to zero (including the absorption) or calculating it from the 
fitted normalisation and temperature of the DISKBB model. The ad- 
vantage of using the second method is that the uncertainties on the 
disc flux are able to be estimated using error propagatiorQ and so 
this method is used for the remainder of this work. 

However, the decay in the disc temperature as the disc flux 
decreases is consistent with the model of Sbisc oc T 4 as shown 
by the li nes in Fig. [12] The square points are from Mill er et all 
j2004alfbl [2006V some of which fall outside the limits of this di- 
agram. These points are consistent with T 4 and the highest temper- 
ature one (0.76 keV) also matches the the set of values we obtain 
(see also Fig. 1 13l>. 

Although we plot both major outbursts here, the decay is clear- 
est in the first one as there is a "clean" decay in flux with only 
a small variation in X-ray colour while in the soft state (Fig. [4j. 
The flux decays slightly when the GX 339-4 reaches the soft state. 

12 A quick investigation into the differences between error propagation and 
a parameter space investigation show that the uncertainties from the error 
propagation are around a factor of two larger than those from the parame- 
ter space investigation. Newer versions of XSPEC will have a Monte Carlo 
method for estimating the uncertainties, which will hopefully allow the util- 
isation of first method (Keith Arnaud, private communication.) 
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Figure 13. The Disc Flux and Disc temperature for observations with an 
X-ra y colour of < 0.2. The x-ax is has been extended to show all the points 
from lMiller etaD i2004blaL l20O6h in white squares. The black solid line is 
the theoretical T 4 relation and the red dashed line is the best fit line to these 
observations T 4,62 . The colour scale shows the date since the first detection 
of the disc in the continuum emission. 



There is a small rise and hardening of the spectrum before the 
GX 339-4 returns to the soft state and the flux and disc tempera- 
ture decay gradually. The entry and exit from the soft state are seen 
as "spurs" off the T 4 relation (Fig. [Til. These may be artifacts of 
the fitting process and result from the limitations of the PCA band- 
width. 

Once the second outburst has reached the soft state, the flux 
rises at a fairly steady X-ray colour and then the spectrum hardens. 
After a gap in the light curve the X-ray colour remains relatively 
steady as the flux decays. Both the entry into and departures from 
the soft state (in the middle and at the end of the outburst) are seen 
as "spurs" off the rest of the points, which scatter neatly around the 
T 4 relation. At face value, the implication of these "spurs" is that 
the disc cools onto the T 4 relation, and at the end of the outburst 
heats up again. This increase in the disc temperature as GX 339- 
4 returns to the hard state is also seen in the summed spectra in 
Section|4] 

All of the points outside of the general scatter around the T 4 
relation are from the observations which have a disc but are in the 
intermediate state. Selecting observations with X-ray colours < 0.2 
and disc fluxes of Soi BC > 5 x 10 _9 ergs _1 retains only those 
which scatter close to T 4 (see Fig. [T3l . Fitting these we obtain a 
best fit of 

O rp4.75i0.23 

with a Spearman Rank correlation coefficient of 0.904 and a 
Kendall Tau of 0.740 with a significance of 23.6a (taking into 
account the errorbars in the disc temperature and the disc flux). 
Although this is still steeper than the theoretical expectation, the 
theoretical expectation still is still a good "by-eye" fit. The range 
of disc temperatures is only just under a factor of 2 and the flux 
a factor of 10, and so the correlation can easily be masked by the 
scatter. 
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Figure 12. For Outburst 2 (LEFT) and Outburst 3 (RIGHT): The X-ray flux against the disc temperature (TOP) and the disc flux against the disc temperature for 
the two main outbursts of GX 339-4 (BOTTOM). The colourscale shows the time since the transition into the soft state - blue is early, red is late. The best fit 
lines are shown by the dashed line (T 9 - 44 for the flux plots, and T 4,75 for the disc flux plots - from the observations with X-ray colour< 0.2). The expected 
correlation of Soisc <x T 4 is shown by the soli d line. The points off the top left of these plots are the observations in the hard state which appear to require 
discs. The square points are measurements fromlMiller et al. (2004b a, 200a). 



6.1 Disc Radii 

If the theoretical relation is applicable, then if the inner radius of the 
accretion disc is constant, the temperature should follow a T 4 law. 
The diagrams show that when in the softest part of the soft state, the 
temperature variation is as expected. Therefore we expect that the 
inner radius of the disc is constant. The inner radius of the accretion 
disc can be calculated from the normalisation of the XSPEC model; 



Rin/Rc 



0.677 



10 kpc 



K 

COS i ' 



(1) 



M/Mq 

where n is the normalisation of the disc model in XSPEC and i = 
30° is the inclinati on of the system. T he distance to GX 339-4 is 
taken to be 8 kpc dZdziarski et al.ll2004l) and the mass as IOMq. 
We do not include any uncertainties from the distance and mass into 
the calculation of the uncertainties of R[ n . The variation in the inner 
radius of the accretion disc during the two outbursts are shown in 
Fig-ElTOP. During the part of the outburst where the temperature 
follows a T 4 law, the inner radius of the disc is constant within 
the uncertainties. The disc behaves as a standard thin accretion disc 
during these parts of the outburst. 

However, the same observations which show a departure from 



T also show a change in the inner radius of the disc. The implica- 
tion of these points is that the disc recedes outwards to something 
like the innermost stable orbit at the beginning of the outburst, and 
falls inwards at the end. This is obviously not a physical explana- 
tion or model for the outburst, and so either the model adopted for 
the emission from the disc is not appropriate or the spectral analysis 
is insufficiently detailed to investigate the disc in the intermediate 

states. 

We follow the analysis of GRS 1915+105 in iBelloni et"ai] 
(2000) and plot simple derivatives of the two measured and inde- 
pendent quantities, T and Ri n (see Fig. 1 14b. It is in principle pos- 
sible to calculate the accretion rate from each spectrum using the 
expressions from a standard thin disc, 



M = 



8TYR'i n a T T 4 
3GM 



(2) 



(e.g. lBellonie TaT] ll997h . Therefore any given Ri D and T implies 
an accretion rate, and lines of constant accretion rate as a fraction 
of the Eddington limited rate where we assume an efficiency of the 
accretion flow, r/, of 10 per cent, are shown in Fig. 1 141 BOTTOM. 
The accretion rate is observed to rise at the beginning of the out- 
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Figure 14. For Outburst 2 (LEFT) and Outburst 3 (RIGHT): TOP: The inner radius of the disc in units of R g through the outburst. BOTTOM: The evolution of 
the temperature of the disc at the inner radius with the inner disc radius. The dotted lines correspond to different constant disc accretion rates. The colour scale 
shows the time since the transition into the soft state - blue is early, red is late. 



burst, reaches a maximum when the state is at the highest soft flux, 
and then falls off. In both the outbursts shown here the peak accre- 
tion rate is approximately 0.03AfEdd, where MEdd = Z/Edd/^c 2 . 
The observations which formed the spurs off the T 4 relation have 
a more rapid change in the accretion rate through t he disc. Com- 
paring these two diagrams with the equivalent one in Bello ni et alJ 
(2000), the slope of the observations in GRS 1915+105 are similar 
to th at of the spurs in GX 339- 4. 

iGierlinski & Pond d2004l) include a colour temperature cor- 
rection factor i n their analysis of disc s in ten black hole binary sys - 
tems, following lshimura & Takahard(ll995l) : lMerloni et al.l l2000). 
The effect of this correction is almost always to harden the spec- 
trum. For the same disc luminosity, adding in the colour correction 
increases the disc temperature. Therefore our results may be consis- 
tent with oc T 4 if the temperature we extract from XSPEC includes a 
colour correction effect. To further investigate whether the relation 
is truly shallower than expect ed more data on the disc in the soft 
state is required. We note that lMcClintock et alJ d2007|) also find a 
dev iation from the expec ted T 4 law in the X-ray Nova H1743-322 
and lTomsick et al.l d2005h present a detailed correlation between L 
and T with deviations from the expected relation for 4U 1630-47. 

Comparing Figs.|7]and[TT] as the disc temperature decreases, 
the equivalent width increases. Excluding the observations of the 



disc in the hard state and fitting the correlation gives 
EW <x T" 3 - 3 , 

with a Spearman Rank correlation coefficient of —0.704 and a 
Kendall Tau of —0.494 with significance of 9.48a. 



7 DISC FRACTION LUMINOSITY DIAGRAM 

In a study whose ai m was to compare AGN and X-ray binaries, 
iKQrding et all (2006) construct a more general version of the HID. 
The location of points in the HID depends on the total luminosity of 
the system and the strength of the non-thermal to thermal emission. 
In binaries, both the power-law and disc components are observed 
in X-rays, whereas in AGN the disc peaks in the UV band. As a 
result the HID for AGN as estimated from X-ray data alone would 
be unlikely to give any insight into the state of the source, and UV 
data from a larg e number of AGN is scarce. 

Therefore Kord ing et ail d2006h generalised the HID by using 
the total luminosity and the disc fraction which is calculated from 
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Figure 15. The Disc-Fraction Luminosity diagram for GX 339-4 for the two 
major outbursts. The colour scale shows the date of the observation relative 
to the start of the outburst - black is early, white is late. The Outburst 2 is 
shown by circles, Outburst 3 by triangles. 



The Disc Fraction characterises the relative strengths of the disc 
and power-law components and also remains finite if either of the 
components approach zero. 

We adapt this characterisation to the spectrum of GX 339-4. 
The model fits throughout the outbursts of GX 339-4 allow the disc 
and power-law fluxes to be calculated, 5di sc from 0.001 — 100 keV 
and Spl from 0.1 - 100 ke\0 This allows us to calculate a Disc 
Fraction Luminosity diagram (DFLD) for an X-ray binary for the 
first time. To calculate the luminosity of the different components 
we take the distance to GX 339-4 to be at 8kpc. The DFLD is 
shown in Fig. \TE\ The colour scale tracks GX 339-4 through the 
diagram during the outburst, and is linked to the time since the 
begin ning of the outburst (similar to Fig.|4j. 

iKording et al simulate a DFLD for 100 objects based 

on previous work on the evolution of the disc and power-law com- 
ponents during an outburst (see their Fig. 10). Reassuringly, even 
for only two outbursts from one object, Fig.[l5]looks similar. 

We note that as we only use X-ray observations in this study, 
this DFLD is limited, as in quiescence a disc is again seen. How- 
ever, this disc is cooler, and so is best observed in the UV and it is 
not seen in the X-ray. This re-appearance of the disc is hinted at by 
the softening of the low flux hard state tail in the HID (Fig. [2}- Fu- 
ture outbursts where the decays are studied in the Ultra-violet will 
reveal the shape of the DFLD for all states that XRBs exhibit. 

The shape of the DFLD is subtly different to the HID. As the 
disc flux and temperature decay, the X-ray colour remains almost 
constant, before very rapidly hardening at the end of the outburst. 
This is a result of the X-ray colour being an indirect measure of the 
disc flux. The Disc Fraction is a more physical measure, and as a 
result the decay of the disc can be easily tracked as the disc fraction 
slowly reduces towards the power-law dominate state. 

The tracks of the two outbursts through the DFLD is very sim- 
ilar to that in the HID. The hardening and softening of the spectrum 
is seen along with the decay in flux. However, the motion is more 
spread out. Whereas in the HID, the range in X-ray colour was 
from 0.03 — 0.3 while the source was in the soft state, the range in 



the disc fraction is from 0.007 — 0.3. This allows a more detailed 
investigation into the variation of the disc during an outburst. 

There is a large gap in the diagram between the power-law 
dominated section and the area where a disc is detected. This is 
the result of the limited spectral coverage of the RXTE PCA. The 
gradual rise of the disc may be taken as the steepening of the lower 
slope in the broken powerlaw. Only once the curve of the black- 
body is dominant does the disc model become the best fitting, by 
which time the disc fraction is already significant. 



8 RADIO/X-RAY CORRELATIONS 

To investigate correlations of the radio flux within th e DFL D the 
we took the com pilations presented in ICorbel et al] feOOOh and 
Gallo et al. (2004). As the radio observations were not coordinated 
with the RXTE observations we match the two data-sets as best 
possible. If there is a radio observation within two days of an X-ray 
observation, then this radio observation is linked to the X-ray one. 
A two day overlap gave a reasonable number of X-ray observations 
which had corresponding radio flux values. 

There was little data from MOSlEE nor from 2.3 and 
1.5 GHz ATCaB observations. As a result we concentrate on 8.5 
and 5 GHz data as well as the spectral index, a (where Studio oc 
u a ). The data from Corbel et al. (2000) cover Outburst 2, whereas 
the radio data from lGallo et alj d2004l) cluster around the peak of 
Outburst 3. 

As can be seen in Fig. [T6j the highest radio fluxes are ob- 
served at the top of the DFLD. There is a sharp drop in radio flux 
as GX 339-4 drops to the bottom of the DFLD and heads back to the 
hard state. There is also a possibility that the radio flux decreases 
as GX 339-4 moves across the top of the DFLD. 

In the right-hand panel of Fig.[T6]the change of the radio spec- 
trum between the hard and soft states is clear. The stalk of the 
DFLD has a spectral index of a ~ 0.2, whereas the soft state has a 
spectral index of a ~ —0.4 — » —0.6. 



9 SUMMARY 

We have performed a comprehensive and consistent investigation 
into the disc and iron lines detected in the X-ray emission from 
GX 339-4. All the 913 public observations in the 11 year RXTE 
archive were reduced, but only 634 had sufficient number of counts 
and flux to be analysed further. Three types of model were fitted 
to each observation; single powerlaw, broken powerlaw and a disc 
+ powerlaw and the best fitting one chosen. The spectra were also 
tested for the presence of a Gaussian line at 6.4 keV. There were 
four outbursts in the data, but we concentrate on the two best sam- 
pled. The relative variation in flux and X-ray colour between these 
two outbursts are remarkably similar. 

A significant iron line was detected in 400 of the 634 observa- 
tions. There is an anti correlation between the flux and the equiv- 
alent width of the iron line for observations of the soft state (X- 
ray colour < 0.22) of EW oc S^lTiokeV • This is steeper than the 
range in slopes found in AGN. In the hard state the data are consis- 
tent with no correlation. As such there is an effect analogous to the 
X-ray Baldwin effect in AGN present in X-ray binaries in the soft 
state. We compare the powerlaw slope, the line equivalent width 



13 The disc flux is the unabsorbed flux, the powerlaw flux is the absorbed 
flux. 
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Figure 16. DFLD with LEFT 8.5 GHz radio flux in mJy and RIGHT a as the colour scale. Points which only have upper limits on the radio flux are shown by 
triangles. 



and X-ray colour of the spectra. The relati on obtained is simila r 
to that seen in a sample of AGN studied by iMattson et alj 12007b . 
Therefore the behaviour of the line is hysteretical. 

Theoretical arguments indicate that the decay in disc temper- 
ature with flux should mimic that for a black body, Sense oc R 2 T 4 . 
We find that the disc flux rather than the source flux needs to be 
used for this relation to be valid. In the soft state the decay in disc 
temperature matches the expected Sbisc oc T 4 . The best fit, how- 
ever, is steeper at T 4 ' 75 resulting from the comparatively large scat- 
ter and small range in disc temperature. This implies that during 
the decay of the disc temperature in the soft state the inner radius 
of the disc is constant. Departures from a constant inner disc radius 
are non-physical, implying that the model for the emission is not 
appropriate or the spectral analysis is insufficiently detailed during 
the the intermediate states. 

Following the method outlined in iKording et all {2006) 
construct a Disc Fraction Luminosity Diagram for GX 339-4. We 
find that the shape qualitatively matches that produced for AGN. 
Linking this with the radio emission from GX 339-4 the change 
in radio spectrum between the presence and absence of the disc is 
clearly visible. The large gap in the DFLD at low disc-high pow- 
erlaw fractions is surmised to result from the lack of spectral reso- 
lution at low energies. Summed spectra from different parts of the 
HID show that this is indeed the case for the low intermediate state. 
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10 APPENDIX: SPECTRAL FITS 



We show in Table [2] the results of fitting the models adopted for 
this analysis to the summed spectra. We only fit the models used 
in this analysis rather than more complex ones which would be 
appropriate for achieving a good fit to this high signal-to-noise data. 
We also show in Fig.[l7]fits to the spectra which did not include a 
line component, as well as the fits used in Fig. [5] but with the line 
normalisation set to zero before producing the plot. 
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Table 2. 
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Table 1. Spectral Fits 



State 


Outburst 


Restrictions' 1 


Exposure 

(ks) 


Ti 


Break Energy 
(keV) 


r 2 


Disc Temp 
(keV) 


Line width 
(keV) 


Line Signif. 


xi 


dof 


1 - Low Hard 6 


3 


F 
D 


5 x icr 10 -» 1.4 x icr 9 

3050 -» 3230 


213.5 


1.563 ±0.001 


65.5 ± 0.8 


2.179 ±0.018 




0.34 ± 0.08 


14.2 


9.608 


(4390.83/457) 


2 - High Hard c 


2 


F 
C 


2 x 10~ 9 -> 7 x 10" 9 
> 0.7 


89.0 


1.720 ±0.002 


39.4 ±0.19 


2.586 ±0.007 




0.75 ± 0.04 


18.60 


5.023 


(2295.32/457) 


3 - Very High Soft 


2 


F 
C 


7.5 x 10~ 9 -> 1.2 x lO -10 
0.09 -> 0.17 


53.4 


2.242 ± 0.023 






0.871 ± 0.002 


0.63 ± 0.06 


7.83 


0.970 


(443.07/457) 


4 - High Soft 


3 


F 
C 


3 x 10 -9 -> 5 X 10 -9 
0.06 -» 0.13 


26.9 


1.939 ±0.061 






0.779 ± 0.004 


0.58 ±0.13 


3.63 


0.989 


(451.96/457) 


5 - Low Soft 


2 


F 
C 


5 x lO -10 -» 8 x lO -10 
0.08 -» 0.13 


10.8 


2.033 ± 0.010 






0.619 ± 0.005 


0.90 ± 0.06 


8.25 


0.946 


(432.32/457) 


6 - Low Intermediate 


2 


D 


2700 -» 2800 


39.6 


2.131 ±0.026 






0.740 ± 0.009 


0.85 ± 0.03 


13.9 


1.059 


(484.08/457) 



C 0.25 -> 0.37 

3.472 ±0.049 9.54 ±0.13 2.104 ±0.030 - 0.29 ± 0.47 2.10 1.325 (605.64/457) 



a Restrictions (F)lux in erg s 1 cm 2 , (D)ate in MJD and X-ray (C)olour. These are for the selections only. The points themselves occur within these, and not necessarily close to the boundaries. Please see Fig. 
4. 

b Low Hard. Adding in a second break improves the fit to x 2 = 2.061; Ti = 1.645, E\ = 10.0 keV, T 2 = 1.438, E 2 = 50.7 keV, T 3 = 2.023. Allowing a disc to the doubly broken powerlaw is a worse fit 
(x 2 = 2.772), however the disc temperature is high at T = 1.82 keV. 

c High Hard. Adding in a second break improves the fit to x 2 = 3.558. Ti = 1.722, E\ = 29.5 keV, T 2 = 2.106, E 2 = 54.6 keV, T 3 = 2.774. Allowing a Disc as well improves the fit to x 2 = 1.332; 
Ti = 1.615, Ei = 28.4 keV, T 2 = 2.121, E 2 = 56.4 keV, T 3 = 2.799 and T — 1.32 keV. Just adding a disc to the doubly broken powerlaw results in a x 2 = 3.384 and a disc temperature of T = 1.22 keV. 
All line components are highly significant using an F-test on the appropriate comparison continuum model. 



